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Abstract 
> 

Quantum billiards provide an excellent forum for the analysis of quantum 
chaos. Toward this end, we consider quantum billiards with time-varying sur- 
faces, which provide an important example of quantum chaos that does not 
require the semiclassical (H — ► 0) or high quantum-number limits. We analyze 
vibrating quantum billiards using the framework of Riemannian geometry. First, 
we derive a theorem detailing necessary conditions for the existence of chaos 

<**] ■ in vibrating quantum billiards on Riemannian manifolds. Numerical observa- 

tions suggest that these conditions are also sufficient. We prove the aforemen- 
tioned theorem in full generality for one degree-of-freedom boundary vibrations 

£> \ and briefly discuss a generalization to billiards with two or more degrees-of- 

vibrations. The requisite conditions are direct consequences of the separability 
of the Helmholtz equation in a given orthogonal coordinate frame, and they 

C^ ■ arise from orthogonality relations satisfied by solutions of the Helmholtz equa- 

tion. We then state and prove a second theorem that provides a general form for 
the coupled ordinary differential equations that describe quantum billiards with 
one degree-of- vibration boundaries. This set of equations may be used to illus- 
trate KAM theory and also provides a simple example of semiquantum chaos. 
Moreover, vibrating quantum billiards may be used as models for quantum-well 
nanostructures, so this study has both theoretical and practical applications. 
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1 Introduction 

The study of quantum billiards encompasses an essential subdiscipline of applied 
dynamics. Within this field, the search for chaotic behavior is one component 
of a large segment of literature concerning quantum chaos. [13, 10, 9] The ra- 
dially vibrating spherical billiard, for example, may be used as a model for 
particle behavior in the nucleus [30] as well as for the quantum dot microde- 
vice component. [21]. Additionally, the vibrating cylindrical quantum billiard 
may be used as a model for the quantum wire, another microdevice. Other 
geometries of vibrating quantum billiards have similar applications. Moreover, 
vibrating quantum billiards may be used to model Fermi accelerators [4, 20], 
which provide a description of cosmic ray acceleration. The study of quantum 
chaos in vibrating quantum billiards is thus important both because it expands 
the mathematical theory of dynamical systems and because it can be applied 
to problems in chemical and mesoscopic physics. 

Quantum billiards have been studied extensively in recent years. [6, 13, 10, 9] 
These systems describe the motion of a point particle undergoing perfectly elas- 
tic collisions in a bounded domain with Dirichlet boundary conditions. Blumcl 
and Esser[6] found quantum chaos in the linear vibrating quantum billiard. Li- 
boff and Porter[19, 24] extended these results to spherical quantum billiards 
with vibrating surfaces and derived necessary conditions for these systems to 
exhibit chaotic behavior. One of the primary goals of this paper is to generalize 
these results to other geometries. 

In the present work, we derive necessary conditions for the existence of 
chaos in vibrating quantum billiards on Ricmannian manifolds. We prove such 
a result in full generality for one degree-of-freedom boundary vibrations (hence- 
forth termed degree-of-vibration (dov)) and also briefly discuss a generalization 
to quantum billiards with two or more dov. In the "vibrating quantum billiard 
problem," the boundaries of the billiard are permitted to vary with time. The 
degree-of-vibration of the billiard describes the number of independent boundary 
components that vary with time. If the boundary of the billiard is stationary, it 
is said to have zero dov. The radially vibrating sphere and the linear vibrating 
billiard each have one dov. The rectangular quantum billiard in which both the 
length and width are time-dependent has two dov. 

The requisite conditions for chaotic behavior in one dov billiards are direct 
consequences of the separability of the Hclmholtz equation [17] in a given orthog- 
onal coordinate frame, and they arise from orthogonality relations satisfied by 
solutions of the Hclmholtz equation. We also state and prove a second theorem 
that gives a general form for the coupled ordinary differential equations that 
describe quantum billiards with one dov. These equations provide an illustra- 
tion of KAM theory, so they are important for both research and expository 
pursuits. 



2 Quantum Billiards with One 
Degree-of- Vibration 

Quantum billiards describe the motion of a point particle of mass too under- 
going perfectly elastic collisions in a domain in a potential V with a boundary 
of mass M 3> Too- (Though mo/M is small, we do not pass to the limit in 
which this ratio vanishes.) With this condition on the mass ratio, we assume 
that the boundary does not recoil from collisions with the point particle con- 
fined therein. Point particles in quantum billiards possess wavefunctions that 
satisfy the Shrodingcr equation and whose time-independent parts satisfy the 
Helmholtz equation. Globally separable quantum billiards with "stationary" 
(i.e., zero dov) boundaries exhibit only integrablc behavior. That is, the motion 
of their associated wavefunctions may only be periodic and quasiperiodic. Two 
types of quantum billiard systems in which this global separability assumption 
is violated are ones with concave boundary components and ones with compos- 
ite geometry. Both of these situations exhibit so-called "quantized chaos" (or 
"quantum chaology").[13, 7] Perhaps the best-known example of a geometrically 
composite quantum billiard is the stadium billiard[15, 3, 22], which consists of 
two semicircles connected by lines to form a 'stadium.' In the present paper, we 
retain the assumption of global separability but permit the boundaries of the 
quantum billiards to vary with time. 

2.1 Necessary Conditions for Chaos 

In Liboff and Portcr[19], it was shown that any fc-superposition state of the 
radially vibrating spherical quantum billiard must include a pair of eigenstatcs 
with rotational symmetry (in other words, with equal orbital (I) and azimuthal 
(to) quantum numbers) in order for the superposition to exhibit chaotic behav- 
ior. One of the goals of the present paper is to prove the following generalization: 

Theorem 1 Let X be an s-dimensional Ricmannian manifold with (Ricmannian) 
metric g. Assume the Helmholtz operator T = V 2 + A 2 is globally separable on 
(X,g), so that one may write the wave- function ip as the superposition 

ip{x) = J2 a n {t)A n (t)ip n {x), (1) 

n 

where x = {x\, ■ ■ ■ ,x s ) is the position vector, n = (m, • • • ,n s ) is a vector of 
quantum numbers, and 

Mx) = f[fj ni) (*j) (2) 

is a product of s "component functions" /• ' (xj). The parameter a„ is a 
normalization constant, and A n (t) is a complex amplitude. If the quantum bil- 
liard of boundary mass M defined on [X, g) experiences one dov in a potential 
V = V{a), where a describes the time-dependent dimension of the boundary, 



then for any fc-term superposition state to manifest chaotic behavior, it is nec- 
essary that there exist a pair among the k states whose s — 1 quantum numbers 
not corresponding to the vibrating dimension are pairwise equal. (That is, 
for some pair of eigenstates with respective quantum numbers (nfe i; • • • , nfe s _ 1 ) 
and (n' k , • • • , n' k ) corresponding to non-vibrating dimensions, one must have 

n k 3 = n 'k- for a11 3 € {!> • • • . s ~ !}•) 



In words, the above theorem states that given a globally separable vibrating 
billiard, a superposition state of a one dov quantum billiard whose geometry is 
described by an s-dimensional orthogonal coordinate system must have a pair 
of eigenstates with (s— 1) equal quantum numbers corresponding to the station- 
ary dimensions of the billiard's boundary in order to exhibit chaotic behavior. 
For a discussion of the separability of the Helmholz operator, see Appendix 
I. Examples of manifolds on which this operator is globally separable include 
well-known ones such as rectangular, cylindrical (polar), and spherical coordi- 
nates and lesser-known ones such as elliptical cylindrical coordinates, parabolic 
cylindrical coordinates, prolate spheroidal coordinates, oblate spheroidal coor- 
dinates, and parabolic coordinates. [23] Note that the preceeding list does not 
exhaust all possible coordinate systems. (There are others in R 2 and K 3 and 
the preeceding examples may be generalized to manifolds in higher dimensions 
for which separability is retained.) Appendix I includes a general procedure for 
determining if the Helmholz equation is separable for a given coordinate system. 

Applying the above theorem to the radially vibrating spherical quantum 
billiard[19], one finds that rotational symmetry between some pair of eigen- 
statesstates in the superposition is required in order for the system to exhibit 
chaotic behavior. That is, the azimuthal and orbital quantum numbers of two 
of the states in the superposition must be equal. The value of the principal 
quantum number n does not affect the existence of chaos. 

The solution of the Schrodinger equation is of the form[25] 

oc 

EiE n t 
A n a„e ~ip n (r). (3) 

Absorbing the resulting time-dependence (in the phase) into the coefficient A n (t) 
yields 

oo 

V>(r, t) = J2 A n {t)a n (t)i> n (r, t). (4) 

In order to examine the present problem, consider a 2-term superposition 
state of the vibrating billiard in (X, g) . The results for a fc-term superposition 
state follow from considering the terms pairwise. The superposition between 
the nth and qth states is given by 

ipnq(x, t) = a n A n (t)ip n (x, t) + a q A q (t)ip q {x, t). (5) 



We substitute this wavefunction into the timc-dcpcndcnt Schrodingcr equation 

ih ^ t ¥ 1= -^ 2 ^ t) > xex > (6) 

where the kinetic energy corresponding to the Hamiltonian of the particle con- 
fined within the billiard is given by 

(7) 
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V, (8) 

where ai,--- ,a s represent the time- varying boundary components, and the 
walls of the quantum billiard arc in a potential V and have momenta Pj with 
corresponding masses Mj. The particle kinetic energy K is the quantum- 
mechanical (fast) component of the Hamiltonian, whereas the remainder of 
the Hamiltonian — representing the potential and kinetic energies of the billiard 
boundary — is the classical (slow) component in this semi-quantum system. We 
use an adiabatic (Born-Oppenheimer) approximation [6] by only considering the 
quantum-mechanical component K of this coupled classical-quantum system as 
the Hamiltonian in the Schrodingcr equation. The Born-Oppenheimer approxi- 
mation is commonly used in mesoscopic physics. In this analysis, we also neglect 
Berry phase. [32] 

For the present configuration, we assume that V does not depend explicitly 
on time. That is, 

V = V(a u --- ,<*.). (9) 

Note that we are applying nonlinear boundary conditions: 

V>(ai (*),-•• ,o 8 (*)) = 0. (10) 

Taking the expectation of both sides of (6) for the state (5) gives the following 
relations: 

-^ 2 ^ q \ - K{\A n \ 2 , \A q \ 2 , oi, • • • , a s ) 

-|p \ = ih[A n A* q + A q A* n + v nn \A n \ 2 + v qq \A q \ 2 + u nq A n A* + v qn A q A* n ] 

(11) 

In a one dov billiard, a(t) = a\(t) is the only time-dependent boundary 
term (with corresponding momentum P(t) = -Pi(t)), so in this case, the kinetic 
energy is written 

K = K(\A 1 \ 2 ,\A 2 \ 2 ,a), (12) 



Ipnq 

8ih.. \ 

ih(tp 



where we use the notation A\ = A n , A2 — A q . The potential energy is given by 

V = V{a). (13) 

In this case, there is a single momentum term in H given by P 2 /2M. Liboff and 
Porter[19] showed for the radially vibrating spherical billiard that if ip n and ip q 
do not have common angular-momentum quantum numbers, then \x nq — /i qn — 
0, where the coupling coefficient fi nn > is defined by 

Vnn' = Mnn'-- (14) 

a 

We show that the vanishing of these coefficients in any one dov quantum bil- 
liard implies non-chaotic behavior of a given superposition state of that billiard. 
Without such cross terms, one observes that Aj is a function of only Aj and a: 

A j = Xj (A j ,a). (15) 

Therefore, \A 3 ;(t)| 2 is a function only of a(t), and so the present system has the 
Hamiltonian 

H(a,P) = K(a) + £-F + V(a), (16) 



2M 



where 



d 
P = -ihV a = a ■ V = -ih— (17) 

da 

is the momentum of the billiard's boundary. The symbol V represents the 
component of the gradient in the direction a. In spherical coordinates, for 
example, we identify a with the unit vector in the radial direction, and thus V a 
represents the component of the gradient in the radial direction. 

The Hamiltonian (16) describes an autonomous single degree-of-frccdom 
system, which corresponds to a 2-dimensional autonomous dynamical system, 
whose non-chaotic properties are well-established. [12, 29] We therefore conclude 
that at least one of the coupling coefficients fi nq or /i qn must be non-zero for a 
quantum billiard with one dov to exhibit chaotic behavior. Wc show below for 
separable systems (see Appendix I) that the condition for the coupling coeffi- 
cients fXnq to vanish is a consequence of the orthogonality of the superposition's 
component functions fj(xj) (see equation (45) in Appendix I). 

In the case of a fc-term superposition, one considers the coupling coefficients 
{/J'nq} °f each pair of cigenstates in the superposition. If any one of these 
coupling terms is non-zero, then one expects the system to exhibit chaotic be- 
havior. Indeed, the fact that the coupling coefficients do not vanish implies that 
one obtains a 5-dimensional dynamical system (which is really a two degree- 
of-frccdom Hamiltonian system in disguise). One observes numerically that no 
matter which two terms one considers, the resulting dynamical system behaves 
chaotically for some set of parameters and initial conditions. Note that the 



above theorem does not hold for two dov quantum billiards, because if one con- 
siders an inseparable potential such as the anharmonic potential, then one has a 
two degree-of-frecdom Hamiltonian system even for cases in which the coupling 
coefficient vanishes. In the next subsection, we discuss the technical details of 
the proof of this theorem. 



2.2 Orthogonality of the Component Functions 

Using the method of separation of variables (again see equation (45) in Ap- 
pendix I) on the Helmholtz equation for a system with s degrees-of-frccdom, one 
obtains s boundary-value problems to solve. (Note that the ordinary differential 
equations for fj are Sturm-Liouville problems.) The solutions to such problems 
may be expressed as eigenfunction expansions [8] . The orthogonality properties 
of the resulting cigcnfunctions are well-known. For each j, 

(#' i /;*)=*»,»;■ as) 

When taking the expectation of the right side of equation (5), the orthogo- 
nality relations satisfied by the s component functions fj play an essential role. 
In the following discussion, fb and mb denote fixed boundaries and movable 
boundaries, respectively. When calculating the expectation, one must integrate 
with respect to all s variables to see when the inner product in (18) is non-zero. 
In particular, this inner product is present in each of the cross terms for the 
(s — 1) fb variables, so those terms vanish unless rij = n'j for each of these 
(s — 1) variables. By the separability of ip, the s-dimensional expectation inte- 
gral is expressible as the product of s one-dimensional integrals, so each term 
includes a prcfactor that consists of the product of (s — 1) inner products. Us- 
ing the Chain Rule, one finds that a variable whose corresponding boundary is 
time-dependent ("m& variables") will manifest differently in the calculation. [19] 
Since the fb variables must have corresponding symmetric fb quantum numbers 
(rij = n'j for all j € {hi, • • • , fc s -i}) for a two-state superposition to have non- 
zero coupling coefficients {/i n <j}, and since we showed above that there must be 
at least one such cross term for a one dov quantum billiard to exhibit chaotic be- 
havior, there must exist a pair of eigenstates whose (s — 1) fb quantum numbers 
are equal. This completes the proof of Theorem 1. 

3 Quantum Billiards with Two or More 
Degrees-of-Vibration 

We now generalize the above results to quantum billiards with vibrations of 
two or more clcgrecs-of-frecdom. Suppose that £ of the s boundary components 
are time-dependent and also suppose that the Hamiltonian is separable: 

« 
H(a u --- ,a^,P u --- ,P i ) = J2H j (a j ,P j ). (19) 



For H to be separable, one requires that both the billiard's potential V(a\, • • • , aj) 
and the kinetic energy i4T(|Ai| 2 , • • • , |Afc| 2 , fli, ■ ■ ■ , a^) be separable in the same 
sense as the Hamiltonian. For some systems, such as the vibrating rectangular 
parallelepiped, the kinetic energy is separable. For others, this need not be 
the case. For example, the spherical billiard has kinetic energy K(r, 9, 4>) — 
Ki(r)K2(0, (j)), which is not equal to K\ + K 2 - 

If, in a given superposition, there are no cross terms in the expectation 
(11'), the £ degree-of-frccdom autonomous Hamiltonian above gives rise to a 2£- 
dimensional autonomous system, which, because of the separability, decouples 
into £ two-dimensional autonomous systems, whose non-chaotic properties are 
well-known. If cither V or K is not separable, the system does not decouple. One 
therefore cannot conclude that the system does not have chaotic behavior even 
in the absense of cross terms. A given system is very likely to behave chaotically 
in this event. In the separable case, then, a fc-term superposition state exhibits 
chaotic behavior when the corresponding (s — £) fb quantum numbers must be 
the same for some pair of eigenstates (that is, the ith fb quantum number in one 
state must be the same as the ith fb quantum number in the other state of the 
pair. Here, i runs over all (s — £) fb quantum numbers). The proof is entirely 
analogous to the one above. 

4 Differential Equations for One 
Degree-of-Vibration Quantum Billiards 

In the present section, we derive a general form for the coupled differential 
equations describing quantum billiards (in a separable coordinate system) with 
one dov and a nonvanishing coupling coefficient (i nq . The resulting system of or- 
dinary differential equations behaves chaotically. Indeed, numerical simulations 
indicate chaotic behavior for some choices of initial conditions and parameters. 

Theorem 2 Consider a one dov quantum billiard with the same geometric con- 
ditions as in Theorem 1. Let the point particle inside the billiard be of mass 
mo, the mass of the billiard's boundary beM> mo, and the surface potential 
of the billiard be V = V(a), where a = a(t). For a two-term superposition, if 
the ith fb quantum number is the same in both states [where i runs over all 
(s — 1) of these numbers], then the system of differential equations describing 
the evolution of the superposition state has the following form in terms of Bloch 
variables x, y, z (defined below) [2], displacement a, and momentum P: 



x = — 



ujqV _ 2[i qq ,Pz 
a 2 Ma ' 



(20) 



(21) 



,-> 



= = T (22 » 



and 



d= M< < 23 > 



p_ dV , 2[e+ + e^(z-p qq ,x)} 
9a a 3 



In the above equations, 



_£,<-£ 



and 



"o = ^-=-^. ( 25 ) 



e ± - ^^, (26) 



where e q and e ? / (g ^ g') are the coefficients in the kinetic energy corresponding 
to the mb quantum numbers. Additionally, x, y, and z represent (dimcnsionless) 
Bloch variables: 

X = Pl2+p21, (27) 

2/ = i(f)2i - P12), (27') 

Z = P22-P11, (27") 

where the density matrix is defined by p qn — A q A* n . [16] 



Before we begin the proof of Theorem 2, note that the differential equations 
describing the evolution of a two-term superposition state are of the above form 
for the linear vibrating billiard [6, 7] as well as for the vibrating spherical billiard 
with both vanishing and non- vanishing angular momentum eigenstates.[19, 24] 
Recall that this system of equations has two constants of motion. They are the 
energy (Hamiltonian) 

H = constant (28) 

and the radius of the Bloch sphere 

x 2 +y 2 +z 2 = l, (29) 

so there are three independent dynamical variables in the set {x, y, z, a, P}. 

We verify Theorem 2 using techniques from Ricmannian geometry It is well- 
known that for s-dimcnsional Ricmannian manifolds, the volume clement dV 
has units of (distance) s , which may include some "prefactors." (For example, in 
cylindrical coordinates, dV = rdrdOdz, where r is the prefactor.) In particular, 
if there are £ "angular variables" (dimcnsionless quantities, like 9 in the above 
example) , there will be a prefactor that includes the term r^ so that the volume 
element has appropriate dimensions. Additionally, in a quantum billiard with 



one dov corresponding to the boundary dimension a(t), the wave- function has 
a normalization factor of order a~ ? , where a corresponds to the number of 
distance dimensions affected by the vibration, which is a different concept from 
the dov. For example, in the radially vibrating sphere, the vibration of the radius 
affects three dimensions, even though this system has one dov. In contrast, for a 
rectangle in which either the length or width (but not both) is time-dependent, 
a single distance dimension is affected, and the dov is also one. When taking the 
expectation of the Schrodingcr equation (6), the normalization prefactor of ip is 
squared, which gives a factor of aT a . We perform s inner products(and hence 
s integrations) in taking this expectation, which yields a factor of a/a in each 
of the cross terms, as was the case for known examples. [11, 1] The diagonal 
terms in the expectation (11) are due only to the kinetic energy, because of 
orthogonality conditions on the wavefunctions ip q and ip q > . 

The evolution equations for A\ and A-i (see (5), etc.) are thus 

2 

iA n = Y J D n jA j , (30) 

where 

{Dnj)= (M ~§)' (3i) 

and \i = n qq i = —[iq'q 7^ is a coupling coefficient (proportional to v qq >) for the 
cross term A q A*, corresponding to equation (11). Transforming these ampli- 
tudes to Bloch variables (27) completes the proof of Theorem 2. The calculation 
is exactly as in the radially vibrating spherical billiard. [19, 24] 

The above equations may be used to illustrate KAM theory. [12, 29, 15] 
Toward that purpose, the number of nonresonant tori that have broken up 
depends on the initial condition of a given integral curve. One may obtain, 
for example, periodic and quasiperiodic orbits (corresponding to closed curves 
in the Poincare map) as in Figures 1 and 2, local ("soft") chaos (in which 
these closed curves become "fuzzy") as in Figure 3, structured global chaos 
(Figure 4), islands of order in a sea of chaos (Figure 5), and finally global chaos 
(Figure 6). The Poincare sections corresponding to the descriptions above for 
the evolution equations of a one dov quantum billiard have initial conditions 
and parameter values given by x(0) = sin(0.957r) w 0.15643446504, y(0) = 0, 
z(0) = cos(0.95tt) « -0.987688340595, ty = 5, a = 1.25, h = 1, ei = 5, 
£2 = 10, and ji = 1.5. Figures 1-6 are plots for the harmonic potential 

V=^(a-a ) 2 . (32) 



5 Phenomenology 

We now discuss the phenomenology of quantum chaos in the present context. 
In the language of Blumel and Reinhardt [7] , vibrating quantum billiards are 
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an example semiquantum chaos, which describes different behavior than the 
so-called "quantized chaos" that is more commonly studied. Quantized chaos 
or "quantum chaology" is the study of the quantum signatures of classically 
chaotic systems, usually in the semiclassical (ft — ► 0) or high quantum-number 
limits. The observed behavior in these studies is not strictly chaotic, but the 
non-intgrability of these systems is nevertheless evident. The fact that their 
classical analogs are genuinely chaotic has a notable effect on the quantum 
dynamics. [18] In the semiquantum chaotic situation, the semiclassical and high 
quantum- number limits are unnecessary and the observed behavior is genuinely 
chaotic. 

In vibrating quantum billiards, one has a classical system (the walls of the 
billiard) coupled to a quantum-mechanical one (the particle enclosed by the bil- 
liard boundary). Considered individually, each of these subsystems is integrable. 
When they arc coupled, however, one observes chaotic behavior in each of them. 
(Physically, the coupling occurs when the particle confined within the billiard 
strikes the vibrating boundary. The motions of the particle and wall thereby 
affect each other.) The classical variables (a, P) exhibit Hamiltonian chaos, 
whereas the quantum subsystem (x,y, z) is truly quantum chaotic. Chaos on 
the Bloch sphere is an example of quantum chaos, because the Bloch variables 
(x, y, z) correspond to the quantum probabilities of the wavefunction. Addition- 
ally, a single normal mode depends on the radius a(t), and so each eigenfunction 
is an example of quantum-mechanical wave chaos for the chaotic configurations 
of the billiard. Because the evolution of the probabilities \Ai\ 2 is chaotic, the 
waevfunction ip in the present configuration is a chaotic combination of chaotic 
normal modes. This is clearly a manifestation of quantum chaos. Finally, we 
note that if we quantized the motion of the billiard's walls, we would obtain 
a higher-dimensional, fully-quantized system that exhibits quantized chaos. [7] 
In particular, the fully quantized version of the present system would require 
passage to the semiclassical limit in order to observe quantum signatures of 
classical chaos. 

6 Conclusions 

We derived necessary conditions for the existence of chaos in one degree-of- 
vibration quantum billiards on Ricmannian manifolds (Theorem 1). In a fc-state 
superposition, there must exist a pair of states whose fb quantum numbers are 
pairwisc equal. The results of this theorem arise from the separability of the 
Schrodinger equation for a given orthogonal coordinate system as well as or- 
thogonality relations satisfied by solutions of the Schrodinger equation. We also 
discussed a generalization of the previous result to vibrating quantum billiards 
with two or more dov. Moreover, we derived a general form (Theorem 2) for 
the coupled equations that describe vibrating quantum billiards with one dov, 
and we used these equations to illustrate KAM theory. 

We showed that the equations of motion (20-24) for a one dov quantum bil- 
liard describe a class of problems that exhibit semiquantum chaotic behavior. [6, 
7, 19, 24] Unlike in quantum chaology, the behavior in question is genuinely 
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chaotic. Additionally, we did not need to pass to the semi-classical (ft — ► 
0) or high quantum-number limits in order to observe such behavior. From 
a more practical standpoint, the radially vibrating spherical billiard may be 
used as a model for particle behavior in the nucleus[30], the 'quantum dot' 
nanostructure[21], and the Fermi accelerating sphere [4]. The vibrating cylin- 
drical billiard may be used as a model of the 'quantum wire' microdevice 
component. [16, 31] At low temperatures, quantum-well nanostructures expe- 
rience vibrations due to zero-point motions, and at high temperatures, they 
vibrate because of natural fluctuations. Additionally, the 'liquid drop' and 
'collective' models of the nucleus include boundary vibrations. [30] The present 
paper thus has both theoretical and practical import because it expands the 
mathematical theory of quantum chaos and has application in nuclear, atomic, 
and mesoscopic physics. 

Appendix I: Separability of the Helmholtz 

Operator 

Consider an s-dimcnsional Ricmannian manifold with metric coefficients 
{ffn,--- ,9ss} defined by 



2 = 1 



'J 



where x% represents the zth rectangular coordinate and Uj represents the distance 
along the jth axis [33]. The Ricmannian metric is then g = Y[ S j=i9jj- For 
notational convenience, one defines hj = Jgjj , so that ^fg = Il 7 -_i hj- ^ n 
cylindrical coordinates in R 3 , for example, x\ = rcos(8), x 2 = rsin(#), and 
x 3 = z, so that one obtains hi — 1, ft 2 = r, and ft 3 = 1. 

We review the following analysis from Riemannian geometry so that the 
proof of Theorem 1 is easier to follow. To express the Helmholtz equation on 
(X, g), one writes the Laplace-Beltrami operator V 2 with respect to the metric 
9- 

V 2 = lyi_f\_A) (34) 

^/9~[ du 3 \9ji du o) 
If the manifold X is three-dimensional, the Laplacian takes the form 
I 

71 
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hih 2 h 3 



A (v__ _jM + A_ (VIjL) + A (\__ A 

dui\guduij du 2 \g22du 2 ) du 3 \g 33 du 3 
d ( h 2 h 3 d \ d ( h 3 h\ d \ d ( h\h 2 d 



du\ \ fti dui I du 2 V h 2 du 2 I du 3 V h 3 du 



(35) 



We now discuss the separability of the Helmholtz equation V 2 ^ + A 2 ?/; = 
0, which is one of our geometrical conditions. To do so, define the Stackel 
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matrix[5, 23] 



S=($ 



13 ii 



(36) 



where $y = $jj(uj), and the {$r,} are specified by the following procedure. 
Define 



C = det(5) = ^2^n M ji- 

3 = 1 

In the proceeding equation, the (j, l)-cofactor Mji is given by 

M^-C-l^'+MetlMO-ll)], 



(37) 



(38) 



where M(j\i) represents the (j, z)-cofactor matrix that one obtains by consider- 
ing the submatrix of S defined by deleting the jth row and the ith column [27] . 
In three dimensions, Mji take the form 



and 



Mn = 



-M 21 



M 31 
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<to 



$12 $13 
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c 



M 



ji 



and 



— = Y[vj(uj), 



3 = 1 

then the solution of the Hclmholtz equation 
separates: 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



i p= nfj( u j)i 

3 = 1 



(45) 
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where fj solves the Sturm-Liouville equation[26] 

££(*l9 +/ '|>*''- ' is{v " 1 * } (46) 

In (46), b\ = A 2 , and all other hi are arbitrary. For a given Stackel ma- 
trix, this prescribes the {rjj} for which the separability conditions hold. It is 
important to note that for a given metric g, the choice of the Stackel matrix is 
not unique and that for some metrics, there is no Stackel matrix that can be 
chosen and hence no way to separate the Helmholtz operator. Note also that 
the time-dependent Schrodinger equation is separable whenever the Helmholtz 
equation is separable. 

As a special case [corresponding to A = in (44)], Laplace's equation is 
separable whenever 



9jj _ M kl 



9kk Mji 



(47) 



and 



— = Mjtflruim), j,k£{l, ■■■,*}. (48) 

9jj 



2 = 1 



Note that the preceeding condition does not completely decribe the separability 
of the Helmholtz operator, because although the Helmholtz equation is separable 
whenever the Laplacian is separable, the converse is not true. 



Appendix II: Galerkin Approximations 

The method used to obtain nonlinear coupled ordinary differential equations 
for the amplitudes A, amounts to applying the Galerkin method [12, 28] to the 
Schrodinger equation, an infinite-dimensional dynamical system. It has been 
used for many years to study nonlinear reaction-diffusion equations that occur 
in fluid mechanics. It can also be used in the study of nonlinear Schrodinger 
equations (NLS). Our treatment of the linear Schrodinger equation with nonlin- 
ear boundary conditions thus parallels established methods for nonlinear partial 
differential equations. Additionally, the Finite Element Method is also based on 
a Galerkin approximation [14] and one can use such methods in inertial manifold 
theory. 

The Galerkin method proceeds as follows. Consider a partial differential 
equation (possibly nonlinear) whose solution is the function i\). Express ip as an 
expansion in some orthonormal set of eigenfunctions ipi(x), i £ /: 

ip(x) = Y^ Ci(x)ipi(x), x £ X, (49) 
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where / is any indexing set and the coefficients Ci(x) are unknown functions of 
some but not all of the independent variables in the vector x, as in the present 
paper. This gives a countably infinite coupled system of nonlinear ordinary 
differential equations for a(x),i G I. (If the partial differential equation is 
linear with linear boundary conditions, then taking an eigenfunction expansion 
gives constant coefficients Ci(x) = Cj.) One then projects the expansion (49) 
onto a finite-dimensional space (by assuming that only a certain finite subset of 
the Ci(x) are non-zero) to obtain a finite system of coupled ordinary differential 
equations. Thus, a two-term superposition state corresponds to a two-term 
Galcrkin projection. If all the dynamical behavior of a system lies on such a 
finite-dimensional projection, then one has found an incrtial manifold of the 
system. [28] 
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Figure Captions 

Periodicity and Quasiperiodicity I in a one dov quantum billiard 
Periodicity and Quasiperiodicity II in a one dov quantum billiard 
Local Chaos in a one dov quantum billiard 
Structured Global Chaos in a one dov quantum billiard 
Islands in a Sea of Chaos in a one dov quantum billiard 
Global Chaos in a one dov quantum billiard 
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Title: Vibrating Billiards 

Date: Mori Aug 2 13:20:46 1999 

a Range = [ 0.802 98013245033106, 2.36754 9668874172 ]; P Range = [ -8.9403973509933756, 9.5364238410596016 ] 

Initial Conditions: f x, y, z, a, P, time )=( 0.1564344 6504, 0, -0.98768834059499999, 1.059602 64 9006622 6, 4.1059602 64 9006628 

Parameters: ( hbar, M, epsl, eps2, mul2, V0, aO )=( 1, 10, 5, 10, 1.5, 5, 1.25 ) 

Num Pts = 5001; Time Step = 0.0050000000000000001 
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